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The thermal entanglement in a two-spin-qutrit system with two spins coupled by exchange inter- 
action under a magnetic held in an arbitrary direction is investigated. Negativity, the measurement 
of entanglement, is calculated. We hnd that for any temperature the evolvement of negativity is 
symmetric with respect to magnetic field. The behavior of negativity is presented for four differ- 
ent cases. The results show that for different temperature, different magnetic field give maximum 
entanglement. Both the parallel and antiparallel magnetic field cases are investigated qualitatively 
(not quantitatively) in detail, we find that the entanglement may be enhanced under an antiparallel 
magnetic field. 

PACS numbers: 03. 65. Ud Entanglement and quantum nonlocality (e. g. EPR paradox, Bells inequalities, 
GHZ states, etc.), 75.10.Jm Quantized spin models, 05. 50. +q Lattice theory and statistics (Ising, Potts, 
etc.), 03. 67. Lx Quantum computation 



I. INTRODUCTION 

It is well known that quantum entanglement |H S Q 
plays a fundamental role in almost all efficient protocols 
of quantum computation (QC) and quantum information 
processing jj, |f| . In one proposal p for physical imple- 
mentation of qubits, a well localized nuclear spin coupled 
with an electron of a donor atom in silicon plays the role 
of a qubit which can be individually initialized, manip- 
ulated and read out by extremely sensitive devices. In 
another 

proposalHSIlGii, 

the spin of an electron in a 
quantum dot plays the role of a qubit. Long decoherence 
time and scalability to more than 100 qubits are two 
of the important virtues of both the schemes. In both 
schemes the effective interaction between the two qubits 
is governed by an isotropic Heisenberg Hamiltonian with 
Zeeman coupling of the individual spins, namely 



H=JS 1 -S 2 +j(S lz + S 2z ). 



(1) 



At extremely low temperatures such a qubit system 
may be assumed to be in its ground state. However a 
real physical system is always at a finite temperature 
and hence in a mixture of disentangled and entangled 
states depending on the temperature. Therefore one is 
naturally led to consider the thermal entanglement of 
such physical systems. The thermal entanglement has 
been extensively studied for various systems including 



isotropic 0> 13 and anisotropic EH Heisenberg 
chains, Ising model in an arbitrarily directed magnetic 
field [lj|, and cavity-QED ll7| since the seminal works 
by Arnesen et aL^H and Nielsen [l^. Based on the 
method developed in the context of quantum informa- 
tion, the relaxation of a quantum system towards the 
thermal equilibrium is investigated [19| and provides us 
an alternative mechanism to model the spin systems of 
the spin-i case for the ap pro aching of the thermal en- 
tangled states [H E H [H [3. It should be noted 
that only the uniform field case is carefully studied in the 
above-mentioned papers. The nonuniform case is rarely 
taken into account. But in any solid state construction of 
qubits, there is always the possibility of inhomogeneous 
Zeeman coupling |20l l2l| . Moreover for perform quan- 
tum computing, it is necessary to control the magnetic 
field at each spin separately 22] . So in the theoretical 
analysis, the nonuniform external magnetic field should 
be included in the model Hamiltonian. Recently, Sun 
[23| and M.Asoudeh [24| investigate the thermal entan- 
glement in the two-qubit spin model with a nonuniform 
magnetic field. But only the spin-^ case is carefully stud- 
ied in the above papers. Zhang et al.|25] only consider 
the uniform magnetic field for spin-1 case. In this pa- 
per, we will investigate the thermal entanglement in the 
two-spin-1 system with a magnetic field in an arbitrary 
direction. Thus we may better understand and make 
use of entanglement in quantum information processing 
through changing the environment. 
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Our paper is arranged as follows: first we will give 
the definition of negativity, the measurement of entan- 
glement. After giving the model Hamiltonian and the so- 
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lutions, we will present our calculation results by several 
figures. Finally, the discussion and conclusion remarks 
will be given. 



II. THE DEFINITION OF NEGATIVITY 

We first introduce the concept of negativity, which 
will be used as the entanglement measure. The Peres- 
Horodecki criterion|26| gives a qualitative way for judg- 
ing if the state is entangled. The quantitative version 
of the criterion was developed by Vidal and Werner j^]. 
They presented a measure of entanglement called nega- 
tivity that can be computed efficiently, and the negativity 
does not increase under local manipulations of the sys- 
tem. The negativity of a state p is defined as 



N(p) = £ \im 



(2) 



where p,i is the negative eigenvalue of p Tl , and T\ denotes 
the partial transpose with respect to the first system. 
The negativity N is related to the trace norm of p Tl 
viaEi 



N(p) 



1 



(3) 



where the trace norm of p Tl is equal to the sum of the 
absolute values of the eigenvalues of p Tl . If N > 0, then 
the two-spin state is entangled. 

The state of a system at thermal equilibrium can be 
described by the density operator p(T) — exp(—f3H)/Z, 
where Z = Tr[exp(— (3H )] is the partition function and 
(i = 1/ksT [ks is Boltzmann's constant being set to be 
unit kg = 1 hereafter for the sake of simplicity and T 
is the temperature). The entanglement in the thermal 
state is called thermal entanglement. 



where J = — |^-, K = -J^j — jfc with t the hopping 
matrix elements, and e = J — K. U s (s = 0,2) repre- 
sents the Hubbard repulsion potential with total spin s, 
a potential U s with s = 1 is not allowed due to the iden- 
tity of the bosons with one orbital state per well, Since 
term e contains no interaction, we can ignore it in the fol- 
lowing discussions and it would not change the thermal 
entanglement. For simplification, J 3> K is assumed and 
the nonlinear couplings is ignored. So the Hamiltonian 
Eq.(4) becomes 



H = J{S X ■ S 2 ), 



(5) 



with an nonuniform external magnetic field in an arbi- 
trary direction, our system is described by 

H = J(S lx S 2x + S ly S 2y )+Bcos[6}S lz +Bsm[9]S 2z . (6) 

in which the neglected exchange coupling term along the 
z-axes is assumed to be much smaller than the coupling 
in the x-y plane. Where S a (a — x,y,z) are the spin 
operator, J is the strength of Heisenberg interaction and 
the magnetic field is assumed to be along the z-axes. 
When the total spin for each site Sj = 1 (j = 1,2), its 
components take the form, 
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III.THE MODEL HAMILTONIAN AND THE 
SOLUTIONS 

The development of laser cooling and trapping pro- 
vides us more ways to control the atoms in traps. In- 
deed, we can manipulate the atom-atom coupling con- 
stants and the atom number in each lattice well with a 
very good accuracy. Our system consists of two wells in 
the optical lattice with one spin-1 atom in each well. The 
lattice may be formed by three orthogonal laser beam, 
and we may use an effective Hamiltonian of the Bosc- 
Hubbard form 28jto describe the system. The atoms in 
the Mott regime make sure that each well contains only 
one atom. For finite but small hopping term t, we can 
expand the Hamiltonian into powers of t and eet|29|. 



In the following calculation we will select J as the energy 
unit and set J = 1. 

To evaluate the thermal entanglement we first of all 
find the eigenvalues and the corresponding eigenstates of 
the Hamiltonian Eq.(6). which are seen to be 



H\ip 3 ) 

mi) 



B+\i; 3 ),H\^t)=-m T \i;t) 



(8) 



H = e + J{S 1 -S 2 ) + K(S 1 -S 2 ) 2 , 



(4) 



where B± = Bcos[6] ± Bsm[6], £ = ^2 + B 2 _, ( = 
'i + B 2 _, m± = \{B+±C)- 

And the corresponding eigenstates are explicitly given 
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by 



IV-i) 

iv>?> 
\*t) 



-1,1)+B_|0,0)-|1,-1)), 



1-1,-1), 

|1,1), 
1 



i 



1,0> + S±|0,-1», 



(|0,1> + S±|1,0», 



^l + R 2 ± + {l±B_R ± y 



- 1, 1) ± R±% 0) + (1 ± B_i2±)|l, -1}]. 



(9) 



r 



where R± = ±5_ + £ and S± = \{B_ ± (). Here \x,y) 
(x = 1,0,-1 and y = 1,0,-1) are the eigenstates of 
Si z S2z- The density operator p can be expressed in terms 
of the eigenstates and the corresponding eigenvalues as 



1 



P 



;exp[-/3£7,] !*,)<*, 



(10) 



z ^ 

where Ei is the eigenvalue of the corresponding eigen 



states and the partition function is Z = 1 + 2 cosh[/?£] + 
4 cosh[i/3C] cosh[^f±] + 2 cosh[/3B+] . 

For our purpose to evaluate the negativity in what fol- 
lowing we need to have a partially transposed density ma- 
trix p Tl of original density matrix p with respect to the 
eigenbase of any one spin particle ( say particle 1) which 
is found in the basis \x, y) (x = 1, 0,-1 and y = 1, 0, — 1) 
as 
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e PB + J 



sinh[i/3C]S_ 



,«±=-^[e-'/3(C± 

±S_ ( 1 - cosh [/3£])-£sinh[/3£] 



where M = cosh[i/3C] 

B+)(-l + 6^)], U± = ^-^-^tJJ-tsmn^tj ^ w± 

£(1 + cosh[/?£](l + S2) ± ^B.sinh^]), Q± = 

l e i/3( C ±S +)(1 + gejjf^W 
2 e ^ C ' 4+S_(S_+C) ' 

We perform the numerical diagonalization of the den- 
sity matrix and the numerical results of the entanglement 
measure N is presented. In Fig.l, we give the contour 
of negativity for different temperature with respect to B 



and 9. From Fig. 1 , we can see that the evolvement of neg- 
ativity is symmetric with respect to magnetic field. The 
maximum negativity arrives at the point B — when 
T = 0.05 and T = 0.2. With the increasing temperature, 
the area of B at which the system can reach maximum 
negativity becomes narrower and even arrives zero for a 
higher temperature (for example T = 0.6 and T = 1.2). 
When the temperature is low, only one peak appears. 
For T = 0.6 and T = 1.2, the double peak structure 




FIG. 1: The contour of negativity for different temperature vs. B and 9. The brighter place means the higher negativity. 



T=0.05 



T=0.2 





T=0.6 



T=1.2 





FIG. 2: The negativity versus B for different temperature, solid curve for 9 = j and dotted curve for 9 = ^p. 



5 



takes place. We can also find that the evolvement of neg- 
ativity is periodic with respect to the polar angle 9 and 
the double peak structure takes place at 9 = (n + |)7r, 
(n = 0, 1, ...). From these figures, we know that the neg- 
ativity gets smaller with the increasing magnetic field 
amplitude. For a higher temperature (say T = 1.2), we 
can see that N arrives zero near B = 0. But N increases 
with the value of B to a peak value for 9 = (n + |)-7T 
cases, after this peak, N will decrease monotonously. 

In order to see clearly the change of the negativity, we 
give the results for 9 = j and 6 = 2f . This corresponds 
to the parallel and antiparallcl magnetic field case. We 
give our calculation results in Fig. 2 (solid line for 9 — f, 
dotted line for 9 = ^), in which the negativity is plotted 
in the whole parameter space at a given temperature, 
and four typical cases are shown. As the temperature is 
low (T = 0.05), we may find that there are two features. 
First, there are three sharp peaks (different from the spin- 
\ case for which only one peak appears, the results for 
spin-i case can be seen from Ref [23]) and the center of 
the middle one locates at B = 0, where the negativity is 
about 1. As we increase the external field B, N rapidly 
decays. When T = 0.6, we can see that the three peaks 
evolve into one. That is to say, the left and the right peak 
disappear, the middle peak gets shorter as we increase the 
temperature. As the temperature is further increased, for 
example T = 1.2, the entanglement is entirely destroyed. 

Compared with the parallel magnetic field case (9 = 
?), when at the low temperature, for antiparallel mag- 
netic field (9 = $2.) case, there is no three-peaks structure 
emerges and only a peak that monotonously decreases 
with the value of B. But the negativity decreases more 
slowly than in the parallel case which means that in the 
strong field region, the parallel field and the antiparallel 
field demonstrate obviously different effects on the en- 
tanglement. We can find that in all parameter space, 
the negativity of the magnetic field with antiparallel di- 
rection is much larger than that of parallel field. These 
results can be seen in these figures. As the temperature 
increase (T = 0.6), the feature of N will be changed. The 
primary peak at T = 0.05 splits into two peaks. For the 
parallel field case, the maximum N appears at B = 0, but 
for antiparallel field case, at the point TV is a minimum 
point. The results is the same with the one which can be 
found in Ref [23]. For T = 0.6, if we apply an antiparal- 
lel field, N will be enhanced more than two times. This 
again illuminates the fact that the well-chosen external 
field can partially weaken the destructive effect of ther- 
mal fluctuation and enhance the entanglement. In other 
words, for a certain temperature, a well-chosen external 
field is helpful for entanglement. At a higher tempera- 
ture (T = 1.2), we can see that N arrives zero at B = 
for the both case. But N increases with the value of B 
to a peak value for the antiparallel case, after this peak, 
N will decrease monotonously. 



IV. CONCLUSIONS 



We investigated qualitatively (not quantitatively) the 
effects of a magnetic field in an arbitrary direction on the 
thermal entanglement in the two-spin- 1 system in terms 
of the measure of entanglement called "negativity" . We 
give results for different temperatures. We find that the 
temperature and the magnetic field can affect the feature 
of the thermal entanglement significantly. At a certain 
temperature, the antiparallel magnetic field is helpful for 
entanglement. In other words, the entanglement may be 
enhanced under an antiparallel magnetic field. 
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